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SECULAR  PERTURBATIONS  If?  THE  ELEMENTS  OF  THE  EARTH'S  ORBIT 

and  astrological  theory  oi~  climate  variations 

By:  Sh.G.  Sharaf  and  N.A.  Budnikova 


Summary 


The  present  paper  deals  with  the  astronomical  theory  of  climate 
variations.  The  tables  of  secular  perturbations  in  the  Earth's  orbital 
elements  for  3x10^  years  time  interval  backward  and  1x10°  years  time 
interval  forward  from  the  epoch  1050.0  are  Riven.  The  secular  course 
of  the  solar  radiation  (in  canonical'  units)  in  the  Earth’s  latitude 
65°  north  and  south  is  illustrated  by  means  of  the  tables  and  diagrams 
for  the  same  interval  of  time  since  1050-0  ar)d  in  the  form  of  the  diagrams 
for  the  past  3x10?  years  and  the  future  lxlOb  years ,  the  "equivalent ' 
latitudes  being  used  instead  of  the  radiation  values  in  canonical  units. 

Hie  problem  of  temperature  cnanges  due  to  the  variations  in  solar  radiation 
arriving  at  the  Earth's  surface  is  discussed,  only  the  influence  of  celes¬ 
tial-mechanic  factors  on  the  Sun's  radiation  being  considered. 


Introduction 


Problems  of  variations  of  the  Earth's  climate  depending  on  changes 
of  astronomical  factors  have  always  been  of  great  interest  to  astronomers 
and  mathematicians.  Mitch,  Kroll,  Boll,  Hargreaves  and  others  have  attempted 
to  determine  the  connection  between  the  climate  variations  and  the  change 
of  the  form  and  position  of  the  Earth's  orbit  in  relation  to  the  Sun. 

However,  only  M.  Milankovitch— a  Yugoslavian  astronomer,  who  devoted  al¬ 
most  his  entire  life  to  thi3  problem  (his  first  studies  in  this  field  ap¬ 
peared  in  1913,  and  the  la3t  in  19>T ) — was  successful  in  creating  an 
orderly  astronomic  theory  on  the  variations  of  the  Earth's  climate  caused 
by  the  solar  irradiation  and  three  factors  of  celestial  mechanics:  inclin¬ 
ation  of  the  ecliptic  to  the  equator,  eccentricity,  and  the  longitude  of 
perihelion  of  the  Earth's  orbit. 

Tbe  solar  constant  and  the  period  of  the  Earth' 3  revolution  around 
the  Sun,  as  values  having  alight  secular  variations,  are  accepted  as  con¬ 
stant  in  this  theory.  In  order  to  characterize  the  climatic  changes  for 
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a  certain  oeriod  of  tine,  M.  Milankovitch  investigated  the  chances  ia  the 
total  solar  radiation  obtained  for  the  sane  period  of  tine  by  a  unit  of 
the  Earth's  surface  at  a  selected  latitude  during  caloric  half-years, 
determined  by  him  aider  the  condition  that  the  quantity  of  heat  obtained 
by  the  uiit  of  area  at  the  latitude  in  any  one  day  of  the  summer  half 
of  the  year  exceeds  the  quantity  of  heat  obtained  by  the  sane  area  in  any 
one  day  of  the  winter  half. 

A  comparison  of  the  amount  of  radiation  obtained  by  the  unit  of  area 
at  the  latitude  if)  during  one  caloric  half-year  of  some  year  in  the  geological 
past  against  the  total  radiation  obtained  by  the  same  area  during  the  same 
half-year  at  the  present,  makes  it  possible  to  Judge  whether  the  given  sur¬ 
face  obtained  more  or  less  heat  during  the  indicated  year  in  the  geological 
past  than  at  the  present;  in  other  words,  to  Judge  about  climate  variations 
for  this  period  of  time. 

M.  Milankovitch  (1931*,  19M)  constructed  tables  and  graphics  of  vari¬ 
ations  in  the  solar  irradiation  of  the  Earth  for  the  600-thousand-year  per¬ 
iod  backward  from  1800;  for  this  purpose  he  utilised  V.  Mishkovich's  (1931) 
computations  of  the  secular  perturbations  in  the  elements  of  the  Earth's 
orbit.  Not  long  ago  Milankovitch 's  results  were  recomputed  'y  Voerkom 
(1958)  with  new  data  on  secular  perturbations  in  the  elements  of  the  Earth's 
orbit  (Brouwer  and  Woerkom,  1950).  He  constructed  the  curves  of  the  summer 
insolation  for  latitude  65°  of  both  hemispheres  for  the  period  of  one 
million  years  backward  from  1950. 

Further  study  of  the  astronomic  theory  of  climate  undertaken  by  the 
Institute  of  Theoretical  Astronorry,  AN  USSR,  on  the  initiative  of  the  All- 
Ubion  Geological  Institute  (VSEGEI),  required  a  computation  of  the  inso¬ 
lation  for  a  longer  period  of  time. 

Preparation  of  the  astronomic  base  of  the  theory  was  the  first  stage 
of  this  investigation.  Without  discussing  in  detail  this  part  of  the 
study,  the  basic  results  of  which  were  published  in  an  article  by  these 
authors  (Sharaf,  Budnikova,  1967),  we  shall  point  out  that  we  reviewed 
and  derived  trigonometric  formulas  of  the  precession  which  consider  the 
second  orders  of  the  eccentricity  and  inclination  of  the  Earth's  orbit. 
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*sd  pracotdiag  fro*  the  values  of  astronomic  constants  in  Brouwer's  and 
Woerkoa's  theory  of  secular  perturbations  in  the  orbit  elements  of  large 
planets  (Brouwer  and  Woerkos,  1^59)  accepted  by  the  International  Astro¬ 
nomical  Uhion,  we  determined  constants  of  Integration  and  coefficients 
in  these  formulas.  We  also  obtained  perturbing  values  for  the  inclin¬ 
ation  of  the.  ecliptic  to  the  equator,  the  eccentricity  and  the  longitude 
of  the  perihelion  of  the  ^rth's  on.it  for  a  39-ad  Xlica-year  period  back¬ 
ward  froa  1959  at  5909  year  intervals. 

Tbe  present  paper  deals  with  the  basic  results  of  the  second  stage 
of  investigations  on  the  astronomic  t  .cor/  of  climte  variations,  i.e.,  with 
the  computation  of  insolation  at  the  u r;er  atmospheric  boundary  for  latitude 
65°  of  both  hemispheres  for  a  39-aillioa-year  interval  backward  froa  1959. 

In  this  ease  we  paid  xuch  attention  to  the  derivation  of  basic  formulas 
of  insolation,  exposure  of  the  cost  influential  periods,  determination  of 
the  auxiliary  constants  and  t.;e  cocparis.cn  cf  our  results  with  the  preceding 
ones.  Wc  should  also  indicate  our  attecot  to  determine  the  connection 
between  the  computed  rate  of  the  radiation  and  the  observed  mean  temperatures 

In  nis  studies  Milankovitch  obtained  theoretically,  after  unavoid¬ 
able  simplifications,  a  simple  relation  tetveen  the  insolation  variations 
at  the  upper  boundary  of  the  atmosphere  and  the  temperature  variations  at 
the  mean  elevation  of  the  land.  According  to  Milankovitch  the  variation  in 
temperature  is  directly  proportional  to  the  variation  in  the  insolation. 

Some  authors,  both  in  the  USSR  and  other  countries,  think  that  the  propor¬ 
tionality  coefficient  obtained  by  Milankovitch  is  strongly  overstated  and 
the  insolation  variations  caused  by  the  factors  of  celestial  mechanics, 
therefore,  are  not  of  great  significance  in  the  history  of  climate.  Usually 
these  authors  refer  to  Simpson's  study  (loliO).  We  were  able  to  prove 
that  Simpson's  study  contains  contradictory  assumptions,  and  his  conclu¬ 
sions  that  the  temperature  variations  do  not  exceed  1.5  -  2°.0  are  not 
quite  true. 

Our  computations  showed  a  good  conformity  between  the  summer  vari¬ 
ations  of  the  temperature  coefficient  and  the  Milankovitch  data.  In  the 
case  of  the  winter  half-year  the  coefficient  proved  to  be  twice  smaller 


-3- 


TC-1681 


we  rvn 


v?-» 


than  for  the  suirrer  half  (Milankovitch ' s  proportionality  coefficient  is 
the  sane  for  both  halves  of  a  year).  On  the  basis  of  this  data  it  is 
possible  to  assert  that  during  the  last  3  million  years  of  geological 
past  the  summer  temperature  could  have  increased  cr  decreased  by  S-6°, 
and  the  winter  temperature  by  2-3°. 

This  paper  concists  of  three  chapters  and  supplements. 

Chapter  I  discusses  theoretical  problems  of  the  astronomic  theory 
of  climate  variations. 

Chapter  II  deals  with  computation  of  the  summer  and  winter  insolation 
for  latitude  65°  of  both  hemispheres  for  a  30-million-year  time  interval 
of  the  geological  past  and  1  million  years  of  the  geological  future.  It 
also  has  some  tables  and  series. 

Chapter  III  is  devoted  to  critical  examination  of  Simpson's  study  and 
to  determination  of  the  relation  between  the  insolation  variation  at  the 
upper  boundary  of  the  atmosphere  and  the  temperature  variations  at  the 
mean  elevation  of  the  land. 

The  supplenents  contain  the  graphs  illustrating  the  secular  rate 
of  insolation. 

The  authors  express  their  deep  gratitude  to  A.V.  Khabakov,  G.S.  Ga- 
neshin,  1. 1.  Krasnov  and  V.A.  Zubakov,  of  the  VSEGEI ,  for  their  attention 
to  this  study,  and  to  R.K.  Sardiyev,  of  the  Institute  of  Theoretical 
Astronomy  A.”  USSR,  for  his  help  in  computations  and  graph  preparations. 
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CHAPTER  I 


BASIC  IDEAS  OF  THE  ASTRONOMIC  THEORY 
OF  CLIMATE  VARIATIONS 

1.  Quantity  of  Solar  Radiation  during  Astronomical  Seasons  of  the 
Year 

If  we  designate  the  mean  amount  of  the  radiation  obtained  by  a  unit 
of  the  earth  area  at  the  latitude  ?  during  the  time  interval  from  t^  to 
t2»  by  W  and  the  mean  annual  rate  of  irradiation  of  a  unit  of  area  at  the 
same  latitude  by  w  ,  then 

t. 

W  =  (  wdl. 

(1) 

M.  Milankovitch  (1939)  gives  the  following  formula  for  determining  w 

w  —  -~r-  sin  ?  sin  3  —  cos  o  cos  3  sin  y0). 


Where  JQ  is  the  solar  constant, 

?  -  latitude  of  the  place, 
p  -  radius-vector  of  the  Sun, 

3  -  Sun's  declination, 

cos00  =  -tg?tg3, 

(2) 

-  and  +  'ro  determine  the  position  of  meridians  which  separate  the  illumi¬ 
nated  part  of  the  parallel  f  from  the  dark  part.  The  Sun.  is  located  at  the 
meridian  o  — 0  .  The  illuminated  part  of  the  parallel  ~  will  be  equal  to 

2Vo  • 

Let  us  re-write  formula  (l)  in  the  form 


W 


1 


sin  3  —  cos  ?  cos  3  sin  y0)  dl. 


•i 


(3) 
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With  the  aid  of  relation 


(3') 


where 

a  =  1  is  the  major  semiaxis  of  the  Earth's  orbit, 

X  -  Sun's  longitude, 

T  -  period  of  the  Earth's  revolution  around  the  Sun, 
e  -  eccentricity  of  Earth's  orbit; 

now  let  us  convert  from  the  independent  variable  t  to  the  variable  >.  in 
formula  (3).  The  new  limits  of  integration,  corresponding  to  t^  and  tp» 
will  be  and  X"  .  Then 


x* 


W 


=  JL  f  — T—i —  (-j(i  sin  -s  sin  ?  —  cos  s  cos  o  sin  y0)  rfX. 

J  \;1  —  c- 


(M 


In  order  to  express  0„  and  0  by  X  ,  we  use  formulas 


sin  D  —  sin  t  sin  X,  cos  %  —  — V  t?f'- 

Here  e  is  the  inclination  of  the  ecliptic  to  the  equator. 

In  the  stated  theory,  values  T  and  JD  are  considered  constant.  The 
eccentricity  e  and  the  inclination  of  the  ecliptic  to  the  equator  e  vary 
within  narrow  limits  (e  from  0  to  0.067,  e  from  22°.068  to  2U°.568)  (Sharaf, 
Budnikova,  1967).  Consequently,  in  the  integration  of  (U)  we  can  disregard 
the  second  power  of  the  eccentricity,  assume  expand  W  in  series 

in  powers  of  and  take  into  consideration  only  the  first-order  terms  with 
respect  to  At  .  Then 


?.* 

> 


(5) 


W0  is  determined  according  to  formula  (U)  in  which  T,  J0,  e»  and  e 
are  considered  as  independent  of  X  in  the  integration,  or,  in  other  words 
independent  of  time. 


-6- 


TC-1681 


Considering  that 


Otu  ,  i)  o  .  .  .. 

— -  =  0,  —  =  c  Lit  s  tir  c, 
rfyo  c»s  °  °  ' 


we  obtain 


dV/  TJ  o 


”  r.  -s'i  —  c- 


1* 


^  (yu  sin  <?  cos  0  —  cos  O  sin  2  sin  O0)  ctg  s0  tg  Zdi.  -f 


L*' 

v 

;  •  \  .* 


^  (C<0  sin  »  sin  o 4-  cos  ?  cos  2  sin  O0)  <2A 


(6) 


Formulas  ( U)  and  (6)  show  that  with  the  replacement  of  o  by  — s 
and  A  by  A-f  ^  (correspondingly  and  are  changed  to  "  and 

*/"-{- k)  the  value  W  does  not  change.  It  follows  hence  that  the  amount 
of  radiation  obtained  by  the  latitude  «?  for  a  period  When  the  Sun  passes 
over  a  segment  of  the  ecliptic's  arc  A"  — A'  is  equal  to  the  amount  of 
radiation  obtained  by  the  latitude -»  during  the  time  interval  when  the 
Sun  passes  over  the  arc  from  A'---  to  A"-^-s 

Let  us  proceed  to  the  determination  of  the  quantity  of  heat  W  obtained 
by  a  unit  of  surface  at  a  given  latitude  -i  during  different  astronomical 
seasons.  Let  Wj,  Wjj,  wjjj»  Wjv  denote  corresponding  quantities  of  ra¬ 
diation  obtained  during  spring,  summer,  autumn  and  winter  by  a  unit  of 
surface  at  latitude  ?  of  the  northern  hemisphere. 

In  determining  values  W  we  should  investigate  two  cases ,  depending 
upon  vhether  the  latitude  ?  is  located  in  the  non-Arctic  or  Artie  zone. 

The  non-Arctic  zone.  In  the  ncn-Arctic  zone  there  is  a  sunrise  and 
sunset  every  day  of  the  year.  During  the  astronomical  spring  the  Sun's 
longitude  varies  from  0  to  90°,  and  during  astronomical  summer,  from 
90°  to  180°.  Consequently,  for  determining  Wj  we  should  replace  A* 
with  0,  A"  with  v  ;  for  determining  Wji  we  should  replace  —A'  with  — 
and  A"  with  -  in  formulas  (U)  and  (6).  In  this  case  it  is  apparent 
that  Wj  =  Wjj.  If  we  designate  the  radiation  for  the  summer  half-year 
by  W8  *  Wj  ♦  Wjj,  we  obtain 
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H\  ---  — —  -----  ^  (v„  sin  o  sin  o  —  cos  o  cos  o  sin  i0)  (!/.  -f- 


n  * 

j  :■  1 

'  !(■/..  ~  i  n  .  ,  .  p  sin i,-, 

rrr.-j  — \  (y0  sin  ~  sin  •>  —  cos  ~  cos  c  sm  o0)  at.  —  cos  o  \ - r1  (2/. 

T.2V1  —  c-’i  j  •  J  cost. 

l.o  0 


(7) 


ClCl-^p,)0. 


For  determining  Wjjj  and  Wjy  the  limits  of  integration  are  replaced 
by  -  and  4 n  for  the  autumn  and  by  and  2-  for  the  winter,  respectively, 
in  the  formulas  (1*)  and  (6).  After  small  transformations,  noting  that 
Wjjj  =  Wjy,  we  can  write  the  following  for  the  radiation  of  the  winter 
astronomic  half-year: 


or 


Ww  =  WII1  + 


w 


1V» 


W„^W. - 


TJt.  r  .  .  , 

'  sm  -o  sm  £a  -L- 

r.  VI—  e-  L 


:su‘ 


sine  cos  s0  (A:)0 


(8) 


The  Arctic  zone.  The  Arctic  zone,  the  latitudes  of  which  satisfy 
the  relations  p  >  -ir  —  e  in  the  northern  hemisphere  ande< — ^-rfin  the 
southern,  has  days  with  sunrise  and  sunset,  without  sunset,  and  without 
sunrise. 

The  Sun's  altitude  in  the  upper  and  lower  culminations  is  determined 
by  the  equalities 


K 

K 


(9) 


Consequently,  during  the  day3  with  sunrise  and  sunset  it  should  be 

/*.><>.  k,<0 
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or 


“(t  -?)<5<-J-r?. 

For  the  days  without  sunset 

*.>0,  5>|  -<f. 

For  the  days  without  sunrise 

/*.<0.  K-tt?- 


Then  in  the  spring  the  declination  o  will  vary:  during  the  -days 
with  sunrise  and  sunset  from  o,  =  0  to  ?2  =  y  —  ?  ;  correspondingly,  dur¬ 
ing  this  time  the  Sun  will  travel  in  longitude  from  /.'  =  (>  to  where 


sin  X 


n 


sin  •j., _ cos  ? 

sin  i  si?:  £  • 


During  the  days  without  sunset  the  declination  6  varies  from  =  —  ~ 

to  B'=e  ;  correspondingly  X'  =  Xj,  X'’  =  y  . 

For  astronomical  summer 


for  the  days  without  sunset  and 


during  the  days  with  the  sunrise  and  sunset.  For  the  astronomical  autumn 


+  ?:  Xf  —  X"  =  «•-{-  X, 
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during  the  days  with  sunrise  and  sunset  and 


0 


1 


X„  X* 


for  the  days  without  sunrise.  For  the  astronomical  winter 


for  the  days  without  sunrise  and 


2 


■  L 


3,=r.O;  ).'  —  2-  —  Xp 


X*  rrr 


for  the  days  with  sunrise  and  sunset. 

By  substituting  the  obtained  values  of  the  limits  X  and  X  in 
equations  (U)  and  (6)  and  considering  that  during  the  days  without  sunset 
—  -r  ,  and  during  the  days  without  sunrise  Vo  — 0  »  w*  obtain  the  fol¬ 
lowing  for  the  latitude  r\  of  the  Arctic  zone  of  the  northern  hemisphere 


IK,  -- — -fr jj  (60 sin «  sin 3 -j-  cos  »eososin*JdX-p«sii!?sin  scosXjj-J- 

rh 

- ...  lr± - ctrr  s  (As)0  \  sin  y  sin  3  -4-  cos  a  cos  o  sin  60)  uX  -}-  -  sin  ?  sin  s  cos  X,  — 

r  K\'i— c*l{iOc  "  [o 

f  cosd,,  n"j 
—  cos  »  l - k-  at-  , 

•  J  COS  S  ’ 

A  11 


W  =  W' - sin  ?  si u  s - sin  ?  cos  * (As)0 

"  *  \'  1  —  c-  5;  V 1  —  c*loO° 


-i 

r/o 


(10) 

(1J) 


As  was  already  indicated,  sinX^g^-  and  A;  is  expressed  in  parts  of  a 
degree.  Let  us  designate  initial  values  of  Ws  and  Ww  by  W°  and  Ww  and 
the  coefficients  of  As  by  AIK,  and  Ali'^  •  Then 

W,*=W*  -j-A!K,As,  \  (12) 

H'.  =  Vn  +  Aiy,As.  | 
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Formulas  (7)  and  (8)  for  the  non-Arctic  zone  yield: 


r*  —  — [  (o  sin  5  sin  l  —  cos  c  cos  o  sin  0, )  d/., 
*  —  C'-  j  ,0  •  ,u/  ’ 


W0  __  IV" - Sin  o  sin  e 

*  *  sVl— C-: 


r  ? .  1 1 

AH'.  =  AH',  -  ilf  cl  5 « (in  -  wy  ) 


For  the  Arctic  zone  we  have  from  formulas  (10)  and  (ll) 


r*i 

II'"  —  — — : —  (  (6,  sin  s  sin  2  --  '-os  »  cos  3  sin  y0)  ■ 

•  *;iVl  —  e-’  J  1 

L" 


■  ^  sin  -  sin  s  cos  /.. 


H/J  —  IV- - 4=^:  sin  s  sin  $ 

*  T.VI-C-- 


AW  -JLrl.rJ  II" _ \  ^  C/Xl  ,  l 

*  «  J 

AH'  —  AH' - -J±=r=-  sin  r  cos 

•  *  lao°  i  1  —  c- 


Aiy„ = nr,  -j-  ~  cta'  £  (U"‘  — 

On  the  basis  of  the  assumptions  regarding  JQ  and  T,  values 
IV0.  M/u,  A1V,  and  AlVv  will  be  constant  for  a  given  latitude.  In  order 
to  determine  them  according  to  formulas  (13)  -  (lo),  we  can  utilize  the 
expansion  of  the  integrands  in  series  in  terms  of  powers  of  small  values , 
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as  this  is  done  by  M.  Milankovitch.  However,  with  an  increase  of  ?  these 
series  converge  slowly;  therefore,  it  is  more  expedient  to  use  one  of 
the  methods  of  quadrature  for  the  integration. 

Formulas  (13)  -  (lo)  were  obtained  for  the  latitude  »  of  the  northern 
hemisphere.  However,  it  was  indicated  above  that  formulas  (1* )  and  (6) — 
from  which  we  proceeded  in  the  derivation  of  formulas  (13)-  (l6) — are  such 
that  with  the  replacement  of  f  by— g  and  ’•  by  X-j--  the  value  W  does  not 
change.  From  this  we  can  draw  a  conclusion  that  a  unit  of  surface  at  any 
latitude  of  the  southern  hemisphere  obtains  as  much  radiation  during  its 
astronomical  summer  half-year  a3  a  unit  of  surface  at  the  same  latitude 
of  the  northern  hemisphere  during  its  astronomical  summer  half-year.  This 
situation  is  also  true  for  the  astronomical  winter  half-year. 

Consequently,  if  we  designate  the  amount  of  radiation  obtained  by  a 
unit  of  surface  at  the  latitude  r  of  the  southern  hemisphere  during  the 
astronomical  summer  and  winter  half-years  by  Ws  and  Ww,  then 


=  (17) 

The  amount  of  radiation  obtained  by  a  unit  of  area  at  the  latitude  ? 
during  the  entire  year  T  will  be 

ivT = ir,  -f  iv  „ = \v]  w*  +  (aiv,  j-  aii'j  hr. 

According  to  formula  (17)  , 


XV T  =  I1V. 

H7*,  IF®,  A!FJt  AH7..  ,  contained  in  formula  (IT),  are  determined  by 
the  same  formulas  (13)  -  (l6).  However,  it  should  be  pointed  out  that 
although  i H and  : Vw  ^  ,  the  thermal  conditions  at  the  same  lati¬ 

tudes  during  the  analogous  astronomical  half-years  in  the  northern  and 
southern  hemispheres  will  be  different.  According  to  Kepler's  second  law, 
the  Sun  passes  over  equal  segments  of  arc  at  different  sections  of  the 
ecliptic  during  different  time  intervals.  Thus,  at  present  the  heat 
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flux  during  a  summer  day  in  the  northern  hemisphere  is  smaller  than  in 
the  southern  hemisphere. 

Let  the  duration  of  the  summer  and  winter  half-years  be  designated 
by  Ts  and  Tw;  the  corresponding  values  for  the  southern  hemisphere  will 
be  Ts  and  Tw. 

Evidently, 


With  an  accuracy  up  to  the  first  order,  with  respect  to  the  eccen¬ 
tricity  of  the  Earth's  orbit  e  we  can  write 


where  II  is  the  longitude  of  perihelion  of  the  Earth's  orbit,  counted 
from  the  point  of  the  vernal  equinox  of  the  date. 

Formulas  (12)  and  ( 17)  together  with  formulas  (13)  -  (l6)  make  it 
possible  to  determine  the  amount  of  radiation  obtained  by  a  unit  of  area 
at  a  given  latitude  of  the  northern  or  southern  hemisphere  during  astro¬ 
nomical  half-years  in  any  year  of  the  geological  past  or  future.  In  this 
case  'I'Y  l-J'Y  H’,„  will  be  functions  of  the  variation  in  the  inclination 
of  the  ecliptic  to  the  equator  e  and  of  values  constant  for  a  given  lati¬ 
tude  and  a  given  half-year. 

However,  a  comparison  of  values  17,,  I7W,  !?,,  I7.r  obtained  according 
to  formulas  (12)  and  (17)  for  a  selected  year  of  the  geological  past 
against  the  magnitudes  of  these  values  for  the  current  period  ’  .  7/j;.,  17",  I?" 
is  not  characteristic  since  these  values  give  the  amount  of  radiation 
obtained  by  the  latitude  ?  during  astronomical  half-years,  and,  as  for¬ 
mulas  (l8)  indicate,  the  winter  and  summer  astronomic  half-years  have  dif¬ 
ferent  duration.  In  addition,  formula  (18)  contains  values  e  and  II 
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which,  in  turn,  change  in  terms  of  time.  Consequently,  the  duration  of 
astronomical  half-year3  also  does  not  remain  constant  in  terms  of  time. 
Thus,  values  !'/„  li'K.  \\\.  !i\  and  »K“.  will  give  the  amount 

of  radiation  for  unequal  intervals  of  timejand  their  variations  in  terms 

of  time  will  not  yield  a  sufficiently  complete  characteristic  of  climate 
for  a  given  latitude. 

For  comparison  we  can  draw  on  other  values ,  namely  the  mean  tension 
of  solar  radiation.  The  mean  tensions  of  solar  radiation  w,  and 
for  the  summer  and  winter  astronomical  half-years  at  a  selected  latitude 
of  the  northern  hemisphere  are  determined  by  the  relations 


and  similarly  for  the  southern  hemisphere 


or  otherwise 


If  we  use  formulas  (12)  and  (18)  and  limit  ourselves  to  the  first- 
order  terms  with  respect  to  the  eccentricity  e  and  variation  in  the  incli¬ 
nation  of  the  ecliptic  to  the  equator  Ae  ,  we  obtain 


*.-T  "i  (‘  -4*»ton)- 


(19) 
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Formulas  (19)  indicate  that  the  variations  in  the  inclination  of 
the  ecliptic  to  the  equator  equally  influence  the  mean  radiation  tension 
of  the  similar  astronomical  half-years  for  a  selected  latitude  of  the 
northern  and  southern  hemispheres.  The  variation  in  the  duration  of 
astronomic  half-years,  the  difference  of  which  is  proportional  to  e  sin 
fl  ,  exerts  an  opposite  influence  in  the  hemispheres. 

Thus,  the  same  latitudes  of  the  northern  and  southern  hemispheres 
obtain  different  mean  radiation  tensions  during  identical  astronomical 
half-years.  The  mean  annual  radiation  tension  will  be 


«*r  =  4-  =  - T  W  + ,|ri  + (AU'‘ 

©r  =  4-  p!';  -f  H'J  4-  (AH\  4-  w  J  *4 


1WU)  A-*]. 


(20) 


2*  Amount  of  Radiation  in  Caloric  Half-Years 

The  variations  of  values  w,  and  — which  depend  on  variations 

of  the  eccentricity,  longitude  of  the  perihelion  and  the  inclination  of 
the  ecliptic  to  the  equator  in  tire — characterize  the  secular  rate  of  the 
irradiation  and  make  it  possible  to  Judge  the  climate  changes  caused  by 
fluctuations  of  the  elements  of  the  earth  orbit. 

In  his  astronomic  theory  of  climate  variations  M.  Milankovitch  prefers 
to  use  the  second  value  for  the  characteristic  of  the  secular  climate 
variations,  namely  the  variations  of  the  amount  of  radiation  obtained 
by  a  unit  of  area  at  a  given  latitude  during  the  caloric  half-years. 

Uhlike  astronomic  half-years,  Milankovitch's  caloric  half-years  are  de¬ 
termined  in  the  following  manner. 

Assuming  that  the  period  of  the  Earth's  revolution  around  the  Sun 
is  a  stable,  constant  value,  i.e.,  the  sidereal  as  well  as  the  tropical 
year  has  an  unchangeable  duration,  Milankovitch  divides  the  year  into 
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two  equal  parts,  ~-Wi  that  one  of  the  half-years,  which  he  calls  summer, 
encompasses  all  the  any s  wnen  the  sum  of  daily  radiation  at  a  given 
latitude  is  larger  than  in  any  of  the  other  (winter)  half-year.  Such 
half-years  are  called  caloric. 

The  difference  between  the  caloric  and  the  astronomic  half-years  con¬ 
sists  in  the  fact  that  during  the  astronomic  half-years  the  Sun  passes 
along  180°  of  the  arc  of  the  ecliptic  and  in  the  summer  astronomic  half- 
year  the  duration  of  any  day  is  longer  than  the  duration  of  any  day  of 
the  winter  half-year.  The  caloric  half-years  last  one-half  of  the  time 
necessary  for  the  Earth  to  revolve  around  the  Sun,  and  during  any  day 
of  the  summer  half-year  the  irradiation  intensity  of  a  unit  of  a  given 
surface  is  larger  than  the  irradiation  intensity  of  the  same  surface 
during  any  day  of  the  winter  half-year. 

The  amount  of  solar  irradiation  for  a  unit  of  area  during  caloric 
half-years  at  a  given  latitude,  determined  for  a  certain  date  of  the 
geological  past  or  future,  can  be  directly  compared  with  the  quantity 
of  solar  irradiation  of  the  same  area,  during  the  same  caloric  half- 
year,  determined  for  the  present,  since  these  insolations  will  refer 
to  the  same  time  intervals. 

In  order  to  determine  the  amount  of  radiation  obtained  by  a  unit  of 
area  at  the  latitude  r  for  the  summer  and  winter  caloric  half-years  (let 
us  designate  these  values  for  the  northern  hemisphere  by  Qs  and  Qw,  and 
for  the  southern  hemisphere  by  Qs  and  Q*) ,  it  is  necessary  to  determine 
the  coordinates  of  the  origin  of  the  caloric  half-years. 

The  total  daily  radiation  IV  ^  ,  obtained  by  a  unit  of  area  at  the 
latitude  -  ,  with  the  Sun  longitude  X  and  the  day  duration  »  will  be 

\VX  (X)  =  ■  (y0  sin  ?  sin  o  cos  ?  cos  sin  y0), 


since  - —  • _ 

1  —  e  cos  (ii  —  X) 

eccentricity 


then  with  an  accuracy  to  the  first  power  of 
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?- 


(21) 


=  1  —  2c  cos  (0  — X). 


Then 

IV,  ().)  =  [1  —  2c  cos  (II  —  >.)1  (-j0  sin  ?  sin  '>  -j-  cos  ?  cos  5  sin  y0).  (  22  ) 

Let  and  *  2  be  the  Sun'5  longitudes  at  the  moments  of  the  begin¬ 
ning  of  the  summer  and  winter  caloric  half-years.  From  the  determination 
of  the  caloric  half-years  at  the  moments  of  their  beginning,  the  diurnal 
totals  of  radiations  should  be  equal,  i.e. 


TV.  ().,)  =  IV,  (X,)7  (23) 

and,  on  the  other  hand,  the  beginning  moments  of  the  caloric  half-years 
should  be  separated  from  each  other  by  a  six-month  time  interval,  equal  to 
We  have 


dt~- 


-  T.  Vi  —  e- 


c?X, 


If  we  substitute  here  the  value  y  from  (21)  and  limit  ourselves  to 
the  first  power  of  eccentricity,  we  can  write 

T“  S  *  =  h  5  1 1  +  2'  "*  (!1  -  ’•>' <r'"  (  2  U ) 

>1 

If  we  solve  equations  (23)  and  (2U)  simultaneously  for  X,  and  X2  we 
can  obtain  the  coordinates  of  the  beginning  of  the  caloric  half-years. 
The  integration  of  (2U)  yields 


n  =  Xs  —  X,  —  ?c  sin  (!J  —  v>)  -;*  '~c  sin  (II  —  X,). 


(25) 
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It  is  apparent  that  *  and  >  differ  from  each  other  by  about  l80°. 
I«t  us  introduce  new  designations:  X,=tX;,  X,=-=r  —  Xr  ,  Then  formula  (25) 
will  be  rewritten  as  follows : 

.  V-f-X'  =  2fIsin(lf  — X')-}-*--iii(n  +  X,)I.  (26) 

Formula  (26)  indicates  that  the  sum  )/_!_)»  will  be  a  small  value, 
on  the  order  of  the  Earth's  eccentricity. 

Let  us  turn  to  formula  (22). 

The  beginnings  of  the  caloric  half-years  can  fall  only  on  days  with 
sunrise  and  sunset. 

Then 


and 


|  cc*s  I  ~  1 1 v*  l<f£|  <  i 


=  7  +  lg  r  lg  2  +  7  lg*  ?  tg1 
sin  6tt  =  1  —  -j  lg*  -»  lg8  Z  —  i-  tg‘ » lg4  Z  —  . 


By  substituting  the  obtained  values  and  sin  in  (22)  we  have 


H\ (X)  -  ~  [1  —  2c  cos  (II  —  >.)) (y  sin  f  sin 2  -i-  y  * *r  «»  ?  cos  S  +  . . 


s  ,  1  sin* ?*ii!*t  ...  n , 

IS  S  COS  0  —  —  •  ■■  ■  i"  ..C  COS  11  ■ 

•  .  ‘  Z  CCS?  CCS  0 

—  xc  sin  II  .sin  s  sin  8  sin  X  —  2c  cos  11  cos  3  cos  S  ccs  X  —  2c  sin  II  cos  ?  cos  5  sin  X^ 


l»',  (X)  =  dji  0.  si  a  o  sin  2  -f  cos  ?  cos  2  -f-  y  •  cos  fl  sin  ?  sin  S  cos  X  - 


(27) 


)/  and  /.'  are  small  values,  as  are  declinations  8'  and  8"  .  If 
we  limit  ourselves  to  the  first-order  values  with  relation  to  c .  X',  X",  S' 

we  obtain  from  (27) 

o 
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Hence 


I  ^(1*) ai fy  siu  j sin i  sin +  (:>»?  —  2ecos  II  cos 
H't  (i.*)  —  »iu  ?  siu  J  *•»  X*  ~r  «®*  r  ~r-c  ca»  1 l cos  r} • 


W.  (X*)  —  IK  (/.')  =  —  f  —  sin  f  sin  i  (.sin  X*  —  sin  >0  -f*  ca»  ^  cas  r*^  —•  & 
**  L*  J 


and,  consequently. 


fin  X*  —  sin  X'  = 


Secos? 
e  tiu  ?  sir.  t 


■  cos  II. 


(28) 


Because  and  X.'  are  small,  their  sines  can  be  replaced  by  arcs. 

Then 


X'  — X'  =  — 


Sc  cos  ? 


r.  sin  f  sin  * 

X#4-X'  sin  11 


cos  n. 


Hence , 


X'«2r(*inn 
>  X'=2c(«in  11 


,  2  cos  o 

I”  r.  sin  c  sin  i 
2  cos  ? 

ft  ^  •  X  •  V  Ntt  * 


cos  11^, 
cos  I  . 


(29) 


(30) 


Values  X'  and  X*  are  functions  of  the  elements  of  the  Earth  orbit 
the  eccentricity,  longitude  of  the  perihelion,  and  inclination  of  the 
ecliptics  to  the  equator,  which  in  turn  are  time  functions.  In  addition, 
X'  and  X"  also  depend  on  the  local  latitude  ~  •  Thus ,  X'  and  X"  vary 
with  tine  and  with  the  change  of  t,h»  latitude  of  the  place.  The  caloric 
half-years  are  not  determined  in  the  equatorial  zone. 

The  arc  through  which  the  Sun  passes  during  a  caloric  half-year  will 
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be  a  segment  of  arc  from  X.  =  X'  to  —  a'  ,  and  during  the  winter 

caloric  half-year  from  l  =■-  «  —  X’  to  X,  — 2*-4-X'  .  In  order  to  determine 
Qg  and  it  is  necessary  to  substitute  these  values  in  the  integration 
limits  in  formula  (U).  We  divide  the  arc  segment  from  X'  to  » —  X’'  into 
three  parts: 

=  Xj  =  U;  ikjesO,  i.?.-:*;  X3  =  s,  a4  -« — V*; 


analogously  the  segment  of  arc  from  r.  —  X'  to  2e-J-X'  is  divided  into  the 
following  parts: 


X|  s-=  t:  —  X',  Xs  =  n;  Xj5»e.  ).,=-^2r;  X.=-2r.  X,  =  2:t-}-X'. 


Then 


k  —  x*)  =  w  (>/, 


0)  + IP.(0.  n)4-ir(n.  s-X") 


or,  since 

IK  (0.  «)r,r/,.  W(r.,  t:  —  X*)  =t  — IV  (r  —  X#,  r),  H'(X',  0)  =  -4V(0,  X').  W  (2e.  2n  --  A1)  —  '.V  Ip,  X'), 


then 


^,*=r/.-ir(0.  X')  -  IK  (n  _  X'.  n). 
ft  X*.  «).f  JK(0,  X*). 


We  designate  W  (0.  A*) -f- \V (r  —  Ap,  r.)r=sK,  then 


Q, 8=3  —  A",  | 

fte=!r„-f  A\  I  (31) 


According  to  formula  (U), 
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K  —  2-7  ,i:;  ?  si»>  5  -r  cos  ?  cos  3  sin  y#)  rfX 


r  (  (V*  sin-  ?  s”»  5  4-  cos  ;  cos  5  sin  y#)  «?>.. 


i-i* 


(32) 


As  was  already  pointed  out  above,  V  and  X"  are  values  on  the  order 
of  the  eccentricity;  then,  as  usual,  by  limiting  ourselves  in  the  inte¬ 
grands  to  the  second-order  values  with  relation  to  and  X*  ,  (32)  may 
be  written 


K  m  !  (r  *!R  ‘ + •»  »)*  -f 


f  sin  ?  sin :  *in  X  4-  cc«  7 j  <?>.. 


(33) 


After  integration  of  (33)  we  obtain 


A  “  {[*” T ?‘:i  ?  ,5n  * **  x  t  ros  ?*J + [— y  sin  ?  sin  1  cos  X  -f  cos  J = 

"*  jJJJ [ — 7  »'•»  ?  sin  *  (ccs  X'  —  cos  X")  4-  n  sin  7  sin  1  v  cos  ?  (X*  -j-  X*)!, 


finally 


by  substituting  here  the  value  X'4-X"  ,  we  have 


K=*JH2g*2.c  s:„  [r 


(3M 


or 


K  rs  M«e  s]||!n. 
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where 


_ 2T/„cos<r 

m  • 


(35) 


Consequently , 


—  me  sin  II,  | 
Qw  =  IV„  -|-  me  sin  II,  ) 

and  for  the  southern  hemisphere 


(36) 


$.  «!•/.+** -in  u.  \  . 

Q0  i=  I i'r  —  me  sin  !  1 .  J 

The  quantity  of  heat  obtained  by  a  unit  of  area  at  the  latitude  ? 
of  the  northern  hemisphere  during  one  caloric  half-year  will  differ 
from  the  quantity  of  heat  obtained  by  a  unit  of  surface  at  the  southern 
latitude  during  the  same  caloric  half-year: 


Q,  —  Q.~ — 2/«csmn,  i 
Qr  —  (fm  —  2 me  kin  II.  / 


The  differences  are  proportional  to  c sin  IT  ;  however,  since  the 
coefficient  a  contains  co3  f  ,  the  differences  will  not  be  the  same  for 
all  latitudes. 


3*  Secular  Variations  in  Solar  Irradiation  of  the  Earth  Surface, 
Caused  by  Variations  in  the  Elements  of  the  Earth's  Orbit 


Formulas  (19)  make  it  possible  to  determine  the  mean  summer  and 
winter  irradiations  obtained  by  a  unit  of  area  at  the  latitude  r  during 
a  time  unit  for  the  moment  t  of  the  geological  past  or  future: 
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«’.=-r,n(>+^i.+-|cSinn). 

®.=4"'i(1  +^SL-'*+T'sii>n). 


(38) 


For  the  initial  moment  t  =  tQ  formulas  (38)  will  be  rewritten 
as  follows : 


—  4<vsin  H0), 

T  V*\>  Uo). 
— —  eo-MaUg). 


i 

1 

(' 


(39) 


The  change  of  value  u-  during  the  period  t  -  tQ  will  be  presented  in 
the  following  form 


*wt  =  y  —  i-lPUc  .sin 111, 

=  4  4-  -i  Wile  Sin  ll] , 

Ate.  =  4  [A!r,Aj  +  4*  xv'?c  sin  n].  ( 140 } 

Alii.  =  4[i»/^s  -  4  Sirt  11]. 

Having  computed  in  a  series  of  moments  the  values  A.v> ,  i;-»  Sui .  '>» 
for  a  selected  latitude  ?  ,  we  will  be  able  to  Judge  the  secular  variations 
of  insolation  at  a  given  latitude.  Since  the  values  l>i\,  l’/t\  IK1' 
contained  in  (Uo)  are  constant  for  a  given  latitude,  tne  variation  of 
isolation  in  terras  of  time  will  depend  on  the  variation  of  values 
li,  Ae.sinlJ  in  terms  of  time. 
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In  order  to  characterize  the  secular  rate  of  solar  irradiation  of 
the  Earth  we  can  also  use  other  relations,  as  was  already  mentioned. 

Formulas  (36)  and  (37)  make  it  possible  to  determine  the  quantity 
of  heat  obtained  by  a  unit  of  area  at  the  latitude  f  during  one  caloric 
half-year  in  the  epoch  t. 


Q,  -- 1  i'-','  -f-  A!  ‘  '.JU  —  inn  sin  !  I. 
(/,  —  IV'l  -f  A! ! 'A:  -f  me  sin  1 1, 


(hi) 


For  the  initial  moment  t  *  tQ  formulas  (hi)  will  be  rewritten  as 
follows 


vK) 

ft  «ir: 


sill  li„. 
sin!!0. 


(1*2) 


The  secular  insolation  variations  during  caloric  half-years  for  the 
time  interval  t  -  tQ  will  be 


A(?,  All  —  rule  sin  1 !, 
A'?,,  —  SV.'wSz  -r  inlr  sin  II, 
—  AH'tAs  -f  mAr  sin  II. 
A(>,r  —  AlK(4.At  —  m.\c  sin  tl. 


(1*3) 


Having  computed,  as  above,  in  a  series  of  moments  A(/s.  -(V  -Q.<-  ^Q., 
we  will  be  able  to  trace  secular  variations  of  the  six-month  irradiation 
for  a  given  latitude.  The  irradiation  variations  in  terms  of  time  will 
depend  on  the  variations  of  A?.  A<*sin!I  in  terms  of  time.  Here,  as 
above,  Ail-,.  All'’,,  and  ra  are  constant  values  for  a  given  latitude  f  . 

It  is  not  necessary  to  discuss  here  in  detail  the  theory  of  secular 
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perturbations  of  the  orbit  elements  of  large  planets. 

In  our  study.  Secular  Variations  in  the  Elements  of  the  Earth's  Orbit 
Influencing  the  Climate  of  the  Geological  Past  (Sharaf,  Budnikova,  1967), 
we  gave  a  short  review  of  the  existing  theories  and  shoved  appropriate  for¬ 
mulas  and  series  for  determining  Is  and  Ac  sin  ft  during  long  intervals  of  time. 
Here,  we  consider  it  useful  to  give  numerical  coefficients  of  series  and 
arguments  for  the  basic  values  which  influence  the  insolation,  namely  for 
the  inclination  of  the  ecliptic  to  the  equator  e  ,  the  eccentricity  of  the 
Earth's  orbit  e,  and  e  sin  0  . 

The  series  for  these  values  will  be  of  the  following  form: 

e  =  h*  -f  ^  /I <  cos  {a.t  -f  b j 
esin  n  =  ^  Ct  sin  !• 

r  =  F„ 4. 2 /’ ,  cos  (/,f J  (^) 


Values  h*,  Ar  t>(  are  given  in  Table  1;  values  cf,  d .  _  in  Table 

2  and  F0.  l\,  g(  ~  in  Table  3.  Tables  1-3  also  contain  values  of  per¬ 
iods  of  corresponding  terms  T^  (in  milleniums). 

As  was  already  pointed  out,  the  insolation  variations  depend  on  vari¬ 
ations  of  c.  csmiH  ,  which  can  be  presented  in  the  form  of  a  sum  of  tri¬ 
gonometric  functions  of  time  with  different  periods  and  coefficients.  Con¬ 
sequently,  insolation  depends  primarily  on  trigonometric  functions  with 
different  periods.  A  greater  influence  will  be  exerted  by  the  terms  with 
greater  coefficients  in  terms  of  absolute  value. 

Thus,  in  the  expression  for  e  there  are  5  terms,  the  coefficients  of 
which  are  larger  than  0°.050  (see  Table  1,  i  =  1  to  5);  the  basic  period 
will  be  a  period  of  about  UlOQO  years.  Because  of  commensurability  of 
these  five  terms,  there  appears  another  period  of  about  200000  years. 

In  the  expression  e  sin  II  attention  should  be  drawn  to  I4  addends  (see 
table  2,  i  s  2  to  5)  the  coefficients  of  which  are  larger  than  0.0050.  A 
period  of  about  20,000  years  will  be  the  most  considerable  period  of  e  sin  n  . 

Six  addends  are  influential  in  the  eccentricity  (see  table  3,  i  = 

10  to  12,  18,  19,  25);  the  coefficients  of  the  remaining  terms  are  smaller 


-25- 


TC-1681 


than  0.0030.  Here  ve  can  speak  about  periods  of  100,000  years,  1*25,000 
years,  and  because  of  commensurability  of  5  term3 — a  1,200  -  1,300  year 
period. 
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CHAPTER  2 


COMPUTATION  OF  SECULAR  VARIATION  IN  THE  SOLAR  IRRADIATION 
OF  THE  EARTH  DURING  A  LONG  PERIOD  OF  TIME 


l».  Tables  of  Secular  Variations  in  the  Elements  of  the  Earth's  Orbit 


We  computed  the  perturbed  values  of  sin  n.  i-co.-FI,  r,  e  for  a  30-million- 
year  period  backward  from  1950  with  an  interval  of  5,000  years.  In  addition 
we  computed  the  perturbed  values  of  A  (e  sin  H)  and  As  for  time  segments 
amounting  to  100,000  years  with  1,000-year  intervals.  Because  of  the 
large  volume  of  tables  the  results  of  these  computations  were  given  in  our 
paper  (Sharaf,  Budnikova,  19&7)  in  the  form  of  graphs,  representing  the  dis¬ 
turbed  values  of  the  eccentricity  of  the  earth  orbit  and  the  inclination  of 
the  ecliptic  to  the  equator.  Here  we  consider  it  useful  to  present  a  port 
of  the  computed  tables,  namely  those  which  contain  the  ~»cular  variations 
of  II  and  e  for  a  period  of  100  thousand  years  backward  from  1950  with  a 
1,000-year  interval  (see  Table  h).  Here 

A  (c  sin  11)  ~  e  sin  if  —  sin  il(>>  ^  —  5  —  V 

The  initial  values  e  sin  n  and  c  for  the  moment  1950.0  will  be 

■  e0  sin  1I0  ==  0.01U ^  23;  '' 4. >7. 

In  Table  5  are  given  A  (r-  si  nil)  and  A=  values  for  3  million  years  back¬ 
ward  from  1950  with  a  5,000-year  interval.  In  Tables  1*  and  5  A^  j-Lir!)  is 
expressed  in  radians,  and  Ac  in  degrees. 
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5.  Auxiliary  Tables, 

In  order  to  detemine  raluee  **'!•  «d  W.  .  formulas  (13) 

and  ( lU )  verm  obtained  in  section  1  for  the  latitudes  of  the  non-Arctic 
tone,  end  (15)  and  (l6),  for  the  Arctic  tone.  The  Simpson  formula  was 
used  in  the  integration.  Values  W*,  H*.  W«.  1W,  and  AH’,  ware  obtained 
for  latitudes  from  30°  to  80°  with  a  2°. 5  interval.  Values  W«,  W*.  H'J  , 
expressed  in  canonical  units  obtained  on  the  assumption  that  the  solar  con¬ 
stant  is  equal  to  1  and  that  time  is  expressed  in  hundred-thousandths  of 
one  year  (T  »  100,000),  are  shown  in  table  6.  Also  given  here  are  magnitudes 
of  values  m  (in  canonical  units)  determined  according  to  formula  (35)  and 
AJF,  and  Air,  which  express  the  variation  of  W#  and  Ww 

with  1°  increment  of  *  .  These  values  were  computed  for  the  intlal 

moment  t0  -  1950.0.  Values  <K*  ft  »  determined  according  to 

formula  (1*2)  in  canonical  units,  are  given  in  Table  T- 

Table  c 
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li 

6.  Computation  of  Insolation  for  a  Thirty-Million-Year  Tine 
Interval  Backward 

Secular  variations  in  the  total  radiation  obtained  by  a  unit  of  area 
at  a  given  latitude  in  a  caloric  half-year,  for  the  time  interval  from  a 
certain  year  in  the  geological  past  till  the  present,  are  determined  ac¬ 
cording  to  the  following  formulas  (see  section  3) : 

&Q'  =  AVI', A*  —  ml  (e  sin  II),  » 

IQ*  =  AU'.fA*  -{-  ml  (c  sin  II),  I 

±Q.  —  AII'.As  -j-  ml  (c  sin  11).  |  (1*5) 

lQm  =  AIV^Ai — mA(e  sin  II).  I 


Values  with  a  dash  refer  to  the  southern  latitude. 
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Table  8  (continued) 
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Table  9  (continued) 


Thuet  haring  all  the  data  contained  in  formula  (1»5),  ve  can  deter- 
alne  the  deviation  of  the  quantity  of  heat  at  a  given  parallel  for  any 
year  of  the  geological  past  from  the  amount  of  heat  falling  on  this 
parallel  at  the  present  moment. 

According  to  formulas  ( U5 ) ,  proceeding  from  Table  6  and  secular  per¬ 
turbations  A  «?  si 't  H )  and  A*  obtained  by  us ,  ve  computed  values  X>., 
and  -(A  for  the  latitude  65®  for  a  30-million-year  period  backward  from 
1950  with  5,000-year  intervals.  Table  8  contains  a  part  of  our  results 
namely  values  aq;,.  A',<,  for  a  3-million-year  period  backward  from 

1950,  with  5,000-year  time  intervals.  Table  9  contains  the  same  values 
for  a  100,000  year  time  period  backward  from  1950  with  the  interval  of 


Wg.  1.  Insolation  during  the  summer  (solid  line)  and  winter  (dotted  line) 
half-years  for  the  latitude  65°  adjusted  to  the  equivalent  latitude 


?  7  *s  ^e  northern  hemisphere,  b  -  is  the  southern  hemisphere.  Time  in 
intervals  of  1.000  years  up  to  1950  is  set  along  the  abscissa  axis. 
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3he  graphs  of  Supplement  1  contain  the  secular  rate  of  the  summer, 
winter  and  annual  insolations  in  canonical  units  for  the  latitude  65°  of 
the  northern  and  southern  hemispheres  for  a  3-million-year  time  interval 
backward. 

According  to  formulas 

Q.*=  ) 

</,-=<r.+  -%■  I  (W) 

X'  1 

ve  can  determine  the  amount  of  irradiation  obtained  by  a  unit  of  area  at 
a  selected  latitude  during  a  caloric  half-year  of  a  certain  year  of  the 
geological  past.  On  the  basis  of  values  SQ,.  -V?„.  obtained  by 

us  and  Table  7»  we  computed  Qs ,  Qv,  Qs ,  Qy  for  the  latitude  65°  encompassing 
a  30-mi Hi on -year  time  period  backward  from  1950,  with  5,000-year  inter¬ 
vals.  We  converted  from  values  Q3  ,  Qy ,  5S ,  5y  to  the  equivalent  latitudes. 
An  equivalent  latitude  is  determined  as  follows.  Let  us  assume  that 
during  a  certain  year  of  the  geologic  past  a  unit  of  area  at  the  lati¬ 
tude  f  obtained  during  one  caloric  half-year  a  total  radiation  equal  to 
Ql*  At  present  the  same  amount  of  radiation  during  the  same  caloric  half- 
y®w  i®  obtained  by  the  latitude  ?,  .  Consequently,  during  the  t]th  year  the 
latitude  *  obtained  as  much  radiation  as  the  latitude  ?,  obtains  now.  The 
thus  obtained  latitude  ?,  is  called  the  equivalent  latitude. 

Fig.  1  shows  solar  irradiation  of  a  unit  of  area  at  the  latitude  65° 
of  both  hemispheres,  reduced  to  the  equivalent  latitude,  for  the  100,000-year 
period  backward  from  1950  (the  computations  were  carried  out  with  an  inter¬ 
val  of  1,000  years). 

The  graphics  of  Supplements  II  and  III  give  the  summer  and  winter  in¬ 
solations  for  the  same  latitude  65°,  reduced  to  the  equivalent  latitude, 
for  a  30-million-year  time  interval  backward  from  1950. 


-1*1- 
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In  addition.  Supplement  II  contains  also  the  curves  of  the  addends 
AU'JU  and  //<A(e  si:i  n ;  of  the  first  and  third  of  formulas  (U5)  and  the 
curves  of  the  perturbed  values  of  the  eccentricity  e.  The  given  curves 
cover  a  3-mi 11 ion -year  period  up  to  1950. 

The  periodicity  which  depends  upon  the  basic  period  of  the  ecliptic 
inclination  to  the  equator  in  1*1  thousand  years  is  clearly  traced  on  these 
curves.  The  peak  amplitudes  follow  the  maximums  and  minimums  of  the  eccen¬ 
tricity.  A  period  on  the  order  of  1200-1300  years  is  also  traced. 

7.  Comparison  of  the  Results  of  Three  Investigations 

In  his  investigations  Milankovitch  (1939;  Milankovitch ,  19l*l)  gives 
a  curve  of  the  summer  insolation  for  65°N  for  a  600 ,000-year  period  up  to 
1850.  In  his  insolation  computations  Milankovitch  was  using  the  compu¬ 
tations  of  V.  Mishkovich  (1931)  who,  proceeding  from  Le  Verrier's  theory 
of  secular  perturbations  of  the  elements  of  orbits  of  large  planets  cor¬ 
rected  for  the  new  values  of  planet  masses,  determined  the  perturbed  val¬ 
ues  of  the  orbit  elements.  In  order  to  determine  the  precession  values 
Mishkovich  used  the  Laplace  formulas  accurate  to  the  first  power  of  the 
eccentricity  and  the  inclination  of  the  Earth's  orbit.  He  used  the  Le 
Verrier-Mishkovich  theory  of  secular  perturbations  in  the  elements  of  the 
Earth's  orbit  and  constants  of  the  Bessel  precession  as  the  basis  for  deter¬ 
mining  the  constants  of  integration  and  coefficients  in  these  formulas. 

In  1950  Brouwer  and  Woerkom  advanced  a  new  theory  of  secular  pertur¬ 
bations  in  the  elements  of  large  planets  (Brouwer  and  Woerkom,  1950).  They 
accepted  new  mass  values  and  took  into  consideration  the  second-order  ef¬ 
fect  caused  by  the  long-period  inequality  in  the  motion  of  Jupiter  and 
Saturn.  Proceeding  from  this  theory,  Woerkom  (1958)  computed  perturbed 
values  of  the  eccentricity  and  the  longitude  of  the  perihelion  of  the 
Earth's  orbit  for  one  million  years  backward  from  1950.  He  obtained  the 
position  of  the  pole  with  relation  to  the  1950  stationary  ecliptic  by 
numerical  integration,  and  then  this  value  was  transformed  for  the  mean 
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position  of  the  non-stationaiy  ecliptic.  In  the  calculation  of  the  pre¬ 
cession  we  took  into  consideration  only  the  mean  value  of  the  position  of 
the  point  of  the  vernal  equinox.  On  the  basis  of  these  results  and  by 
utilizing  the  ,  *s ,  .'-05,  o  values  obtained  by  Milankovitch ,  Woerkom 
(1938)  constructed  the  curves  of  the  summer  insolation  for  the  latitude 
65°  of  both  hemispheres  for  one  million  years  backward  from  1950. 

As  mentioned  above ,  we  obtained  the  insolation  curves  by  proceeding 
from  the  perturbed  values  of  inclination  of  the  ecliptic  to  the  equator, 
the  eccentricity  and  the  longitude  of  the  perihelion  of  the  Earth's  orbit 
obtained  by  us. 

Me  based  the  computation  of  these  values  on  Brouwer's  and  Woerkom's  theory 
of  secular  perturbations  in  the  elements  of  orbits  of  large  planets;  in 
order  to  determine  the  precession  values  we  used  our  trigonometric  pre¬ 
cession  formulas  (Charaf,  buinikova,  106?),  in  which  we  accepted  the 
recent  values  of  astronomic  constants  recommended  by  the  International 
Astronomical  Union. 

In  Fig.  2  are  given  curves  of  the  solar  insolation  for  65°N,  obtained 
by  Milankovitch,  Woerkom  and  us.  As  the  curves  3hcw,  despite  some  dif¬ 
ferences  in  the  initial  data,  the  results  of  the  three  investigations  agree 
rather  well. 


8.  Computation  of  Insolation  for  One  Million  Years  Forward 


The  determination  of  the  insolation  variation  under  the  influence  of 
the  celestial-mechanical  factors  and  for  the  geological  future  is  also  of 
great  interest. 

We  computed  the  perturbed  values  of  the  inclination  of  the  ecliptic 
to  the  equator,  the  eccentricity  and  longitude  of  the  perihelion  for  one 
million  years  forward  from  1950  with  an  interval  of  5000  years.  The  graphs 
of  Fig.  3  and  1*  show  perturbed  values  of  the  eccentricity  and  the  inclin¬ 
ation  of  the  ecliptic  to  the  equator.  Table  10  contains  secular  variations 
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•  sin  H  and  e  for  a  period  of  one  million  years  forward  from  1950,  with 
and  interval  of  5000  years. 

The  values  **(?.•  ~Q .  for  the  latitude  6^°  for  one  million 

years  forwaid  are  given  in  Table  11. 


Table  10 
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Supplement  IV  gives  the  summer  and  winter  insolations  for  the  lati¬ 
tude  65°  of  the  hemispheres,  reduced  to  the  equivalent  latitude,  for  a 
period  of  one  million  years  forward  from  1950. 

Fempl  (1958)  computed  insolation  variations  (AO  ,  AO,,  AO.) 
for  various  latitudes  for  a  period  of  100,000  years  forward.  His  results 
agree  well  with  ours  for  this  period  of  time  and  the  latitude  65°. 


Fig.  1».  Perturbed  values  of  inclination  of  the  ecliptic  to  the 
equator. 

Time  in  milleniums  from  1950  is  given  along  the  abscissa  axis  . 


CHAPTER  III 

RELATION  BETWEEN  THE  SOLAR  IRRADIATION 
OF  THE  EARTH  AND  TEMPERATURE 


In  his  monographs  M.  Milankovitch  (1939,  19bl)  studied  in  detail  the 
influence  of  the  atmosphere  on  the  quantity  of  heat  obtained  by  the  earth 
surface  from  the  Sun  and  derived  the  relationship  between  the  temperature 
at  the  level  of  the  mean  eleVation  of  the  continent  and  the  insolation. 

Both  in  the  USSR  and  other  countries ,  on  opinion  is  expressed  in 
numerous  studies,  devoted  to  the  problems  of  climate  variations,  that 
the  temperature  variations  caused  by  the  insolation  variations  were  con- 
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isiderably  overstated  by  M.  Milankovitch ,  and  the  insolation  variations 
caused  by  variations  in  the  ele:nents  of  the  orbit  and  the  axis  of  the 

rotation  of  the  Earth  are  not  of  great  significance  in  the  paleoclima- 

• 

talogy. 

These  remarks  concerning  the  astronomic  theory  of  climate  variations 
ere  primarily  based  on  a  part  of  Simpson's  study  (19^0)  which  criticizes 
Milankovitch's  conclusions.  Therefore,  we  consider  it  useful  to  discuss 
Simpson's  study  in  detail  and  attempt  to  find  a  connection  between  the 
insolation  variations  and  the  variations  in  real  temperature. 


9.  Variation  of  Temperature  of  the  Earth's  Surface  Caused  by 
Variation  of  Solar  Irradiation  (According  to  Milankovitch) 


M.  Milankovitch  (1939)  assumes  that:  1.  The  surface  of  the  Earth  is 
uniform  (continuous  land)  and  entirely  smooth;  2.  The  Earth's  atmosphere 
is  stationary,  the  air  is  viewed  as  an  ideal  gas;  3.  The  conditions  of 
absorption  and  reflection,  atmospheric  conditions  of  seasonal  variations 
do  not  endure. 

Uhder  these  conditions  and  after  unavoidable  simplifications  Milanko¬ 
vitch  derived  a  fo-mula  yielding  the  relationship  between  the  temperature 
at  the  level  of  the  mean  elevation  of  the  land  and  solar  irradiation: 

c0<  =  l(i  _,i)(]  -f/...;/)ir.  (U7) 

Here:  P  is  the  mean  absolute  temperature  of  the  lower  layer  of  the  air, 
corresponding  to  the  half-year  under  investigation,  at  the  latitude  ? 
during  a  selected  year  of  the  geologic  past  t;  W  is  the  quantity  of  ra¬ 
diation  (in  calories)  obtained  by  a  unit  of  surface  (1  cm?)  at  the  atmos¬ 
pheric  boundary  in  a  unit  of  time  (l  min)  during  the  same  half-year  at  the 
latitude  r  ,  in  the  epoch  t;  0  is  the  Stefan-Boltzmann  constant;  A  is  the 
albedo  (with  consideration  of  the  reflection  of  the  atmosphere);  k  is  the 
coefficient  of  absorption  for  the  sun's  thermal  radiation  in  the  atmosphere, 
M  is  the  mass  of  air  located  above  a  surface  unit. 
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For  the  initial  moment  tQ  the  values  contained  in  formula  ( UT )  are 
designated  by  0.,  ,  Aq,  kc.  Mq,  Wq;  then 

C0j=l(i-,iu)(i-f/.l);./u)nv  ( i*8) 

Assuming  that  M  does  not  change  in  terms  of  time;  that  the  variation 
of  the  k  value  can  be  caused  only  by  an  increase  or  decrease  of  the  quan¬ 
tity  of  vapor  or  carbon  dioxide  in  the  air  and  therefore  k  can  assume  the 
mean  value  k  *  kQ;  and  then  passing  from  the  absolute  temperature  to 
the  temperature  in  degrees  Celsius 

0  =  273°  H-u, 


we  obtain  for  the  difference  of  temperatures  between  the  epochs  t  and  tQ 


Au  = 


i  4-  kryn. 
Si  (273)i 


f(l  _  A)  ]V  —  (1  —  A0)  TFJ. 


(*•9) 


If  we  limit  ourselves  to  the  investigations  of  the  areas  of  the  earth 
surface  which  were  not  buried  under  the  continental  ice  during  the  geological 
past,  we  can  assume  A  *  Aq,  then 


Au  =  /,11V, 


(50) 


here 


1  h,M n  , ,  j  . 


Bgr  assigning  to  the  quantities  contained  in  n  their  mean  values 

yi.~0.40,  /,•„  =  0.0C25,  ^7.  =  3 033.3,  s  =  0.7G .  10*;\ 


we  obtain  Au.=  173.8131' 


.  or  considering  that 


where  is  expressed  in  canonical  units ,  we  have 


Aur^i(A-»t~  A<0- 


(51) 


10.  Connection  Between  the  Temperature  Variation  and  Insolation 

(According  to  Simpson j 


Simpson  (l^O)  assumes  that  the  coefficient  in  Milankovitch's 

formulas  (51)  is  incorrect,  since  it  was  computed  proceeding  from  the 
mean  values  Aq,  kQ,  Mq,  o  ,  i.e.,  values  which  thus  far  cannot  be  deter¬ 
mined  with  sufficient  accuracy.  According  to  Simpson,  because  of  the  com¬ 
plexity  of  meteorological  causes  of  the  climate  variations  only  empiric 
methods  can  assist  in  the  solution  of  this  problem. 

He  suggested  another  method  for  determining  the  connection  between 
the  temperature  and  the  insolation  variations.  Here  is  a  brief  explanation 
of  his  conclusions. 

Let  u^,  Uj  be  the  mean  air  temperatures  of  the  warmest  and  coldest 
months ,  respectively,  at  the  latitude  ?  for  any  year  of  the  geologic  past 
t;  Qg,  Qw  -  total  radiation  (in  canonical  units)  obtained  by  a  unit  of 
area  at  the  latitude  v  during  caloric  half-years  in  the  epoch  t. 

Values  pertaining  to  the  present  time  tQ  are  provided  with  a  super¬ 
script  zero. 

We  de3i<piate 

Y  =  U,  -  H„. 

Ai/j  =  ft]  —  uj, 

Au.,  —  u,  —  !f  j. 

A/f  — :  ;;  —  Jl". 

AV  =  Y  — 
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Simpson  assumes  that  the  mean  temperature  of  the  warmest  month 
can  be  identified  with  the  mean  summer  temperature,  and  the  mean  tempera¬ 
ture  of  coldest  month  can  be  identified  with  the  mean  winter  temperature; 
that  the  difference  between  the  summer  and  winter  insolation  for  any  place 
and  for  any  date  varies  in  the  same  proportion  as  the  difference  between 
the  summer  and  the  winter  temperature,  and  the  ratio  of  variation  between 
the  summer  and  the  winter  temperatures  is  directly  proportional  to  the 
ratio  of  variation  between  the  summer  and  the  winter  variations.  That  is, 


l‘n 


«»  "  II-.,  f 

(53) 

j;  __  v 

v~‘  • 

(5*0 

(55) 

On  these  assumptions  Simpson  obtains  the  following  formulas,  connecting 

the  temperature  variation  with  the  insolation  variation: 

y.»  v 

■'<?'■)  (56) 

By  substituting  here  values  Q'\  (*',  ~Q.  from  Milankovitch's  (1939) 

tables,  and  by  replacing  us  and  Uy  with  observed  mean  temperatures  of  the 
wannest  and  coldest  months  for  the  corresponding  latitudes,  Simpson  com¬ 
puted  and  Af/_.  for  a  number  of  years  of  the  geologic  past  and  con- 

pared  these  values  with  d/r  and  computed  according  to  Milankovitch's 

formulas  (51).  The  comparison  shows  that  and  obtained  according 

to  Simpson's  method  are  about  3  times  smaller  than  the  values  computed  ac¬ 
cording  to  Milankovitch's  formulas.  In  Simpson's  computations  -w/  and 
lu,  do  not  exceed  1.5-2°.0  for  the  deepest  insolation  peaks.  Hence 
Simpson  arrives  at  a  conclusion  that  Milankovitch's  coefficient  p--  is 
highly  overstated. 


-51- 


TC-1681 


However,  this  is  not  auite  so. 

Formulas  (5M  and  (55)  as  well  as  (56)  are  correct  only  when  the 
conditions  discussed  in  section  9  are  observed.  Consequently,  these 
relations  will  be  true  only  with  substitution  of  "fictitious"  tempera¬ 
ture  values  obtained  theoretically.  Naturally,  the  truly  observed  mean 
temperatures  do  not  satisfy  these  conditions.  It'  is  sufficient  to  say 
that  meteorological  conditions  are  different  during  the  summer  and  winter; 
the  snow  cover  changes  the  albedo  value  sharply.  Thus,  it  can  be  ascer¬ 
tained  that  Simpson's  reasonings  contain  not  quite  acceptable  assumptions. 


11 «  Relation  Between  Solar  Irradiation  and  the  Observed  Temperature 


A  number  of  authors  give  tables  of  mean  latitudinal  temperatures. 
Usually  the  mean  temperatures  of  the  warmest  and  the  coldest  months  of 
the  year  are  determined.  With  the  aid  of  Milankovitch 's  formulas  (139, 
p.  171)  it  is  possible  to  pass  from  the  mean  temperatures  of  the  warmest 
and  the  coldest  months  of  the  year  to  the  mean  temperatures  of  the  sum¬ 
mer  and  winter  half-years. 

The  Milankovitch  formulas  yield 

n,  -  0.822;, ^-U.iTiv,  1 

iV.--U.822;;.  —  0.178a..  J  ^7) 

In  T.F.  Betlyoyevay  article  (i960)  a  table  is  given  for  the  monthly 
distribution  of  mean  temperatures  for  the  latitudes  every  5°  from  -70°  to 
+80°.  d.I.  Stekhnovskiy  (1962)  gives  mean  temperatures  of  lands  and  oceans 
for  the  same  latitudes,  also  every  5°. 

On  the  basis  of  T.F.  Batlyayeva's  and  D.I.  Stekhnovskiy' 8  data  we  com¬ 
puted  the  mean  values  of  summer  and  winter  temperatures  for  all  latitudes 
of  continents  and  oceans  (see  Table  12). 

In  Fig.  9  are  given  the  curves  of  distribution  of  the  mean  temperature 
by  tones  for  all  the  Intitules  and  mean  temperatures  of  lands  and  oceans 
at  a  given  latitude. 
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k. (»  i  >  *i  !••••  it.ro  HJtt  I'm*/ 

11.7  i..t  ht.f  j  H  •«  I  t  l>t  «i  l«i 

l. 7  ?.!»  t.#  /  :vi  i  ■*..>»  t,<«» 


In  section  5  (Table  7)  we  determined  the  totals  of  radiation  Qs  and 
Qw»  obtained  by  the  latitude  ■?  during  one  caloric  half-year  for  t  =  t0  = 
1950.0,  expressed  in  canonical  units.  Qs  and  Qy  differ  from  the  mean 
quantities  of  radiations  Ws  and  Wv  only  by  a  constant  multiplier. 

Let  us  assume  that  the  variation  of  the  mean  observed  temperature 
by  latitudes  corresponds  to  the  insolation  by  latitudes,  and  the  deviations 
in  mean  temperatures,  which  depend  on  geographic  conditions,  are  of  an 
accidental  nature.  We  assume  that  such  a  relation  is  also  preserved  for 
the  geological  past. 

Then  the  mean  temperature  and  insolation  can  be  considered  as  con¬ 
nected  by  a  functional  relation  of  the  following  type 

L  ±Q  j  (58) 

here 

iC',  ==<?«-<?“  . 

Qg,  Qv  are  insolations  at  our  selected  latitude  'P.,  = — c  ,  a, ,  bf.  bw 
are  the  coefficients  subject  to  the  determination. 

We  compile  two  systems  of  1?  conditional  equations,  corresponding  to 
latitudes  from  25  to  80°  every  5°,  for  the  northern  hemisphere.  The  other 
system  pertains  to  the  winter  half-year. 

Similarly,  for  the  southern  hemisphere  we  also  compiled  two  systems, 
but  consisting  of  10  conditional  equations ,  corresponding  to  latitudes 
from  25°  to  70°  with  5°  intervals. 

We  substitute  3  variants  of  the  mean  seasonal  temperatures  (of  all 
the  lands  and  oceans)  in  the  right  parts  of  these  conditional  equations. 

We  designate:  as  the  mean  summer  and  winter  temperatures  of 

the  entire  latitude  ;  os  the  mean  summer  and  winter  land  tem¬ 
peratures  at  the  latitude  ;  ,  the  mean  summer  and  winter 

temperatures  of  the  oceans  at  the  latitude  Values  with  a  dash  pertain 
to  the  southern  hemisphere. 
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Fig.  5.  Distribution  of  the  mean  temperature  for  the  entire  latitude 
u,  land  u'  and  oceans  u"  by  latitudinal  zones. 

Tables  13-16  contain  the  conditional  equations  and  the  systems  of 
normal  equations  that  correspond  to  them  (v,  v' ,  v"  are  the  values  of 
discrepancies ) . 

The  solution  will  be  as  follows: 
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*  This  conditional  equation  will  have  only  two  variant#  •  Thli  conditional  equation,  like  the  previous  one,  will 

of  the  right  hand  part.  At  the  latitude  60°S  there  la  have  only  two  variants  of  the  right  hand  part, 

no  land.  Therefore*  we  have  here  two  systems  of  normal 

equations . 


u 


—  2^'i  4-Or007411.ua 
u'  ~  AM  -f  UMOo.MM't 
u\  —  2.3S  4-  o.oococoia 


IvviS,  ) 
i-73,  < 
2.t>5,  j 


(59) 


u.  —  -i5?-'0 Or003943X'«  ±  0?S7,  | 
«;*=—  ii).r*s -f  0.004 1 s-.n.o,  -  i .07 ,  ; 
«'  =  —  8.5S  4-  0.0C‘3:iJi0i(^  ~  4.37,  j 


(60) 


u,  **  —3:07  4-  0:007623.1#,  ±  -O'*,  j 
»  — 4.03  -f  0.006010^,  ~  1 .43.  j 
ti; «  — 3.oo  -f  0.0070$:*^  ±  \  .03,  j 

»  —  Orl"  4-  0:0034303#.  ±  3:4t,  j 
—  —11.80  4-  0.0O2GGO4#,  *  4.21 ,  j 
S.tM  +  0.003200^*2.24.  j 

r 

For  the  purpose  of 'control,  6  additional  systems  of  7  conditional 
equations  corresponding  to  the  latitudes  fron  20°  to  80°  with  10°  inter¬ 
vals  were  compiled;  they  are  of  type 

at  *j*  bt^Q(  —  «4» 
c„  H*  W*> 

c  j  — yt :;s 

In  their  right  hand  parts  we  substituted  the  mean  temperatures  of  the 
entire  latitude,  obtained  from  the  Meynardus  tables  (Simpson,  1Q1*0). 
Following  is  the  solution  of  these  problems: 

u t  =  2?42  -f  OrOoTOi  u2(J.  ±  1  ?o0.  i 

u„*~=~17.3  +0/Oi0CG^::.  J-A  j  (63) 

u f  ■—  —  8.24  -f  0.0032G2A#..  j;  1.04,  j 


(61) 


(62) 
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u 


'  —  —  o.OG  -i-  O.OOjtfOfiX'.  ±  2.21 , 

ug r_--  - 1 7.0  0.0072f;$Xu 5/iS, 

,7  =  _1 1 .9  -i-  0.G0.7  jSuA/V..  3.40. 
2  '  v  . 


(6U) 


The  great  significance  of  the  mean  square  error  for  equations  of  the 
northern  hemisphere  in  the  winter  half-year,  in  the  case  where  the  ocean 
temperature  is  under  investigation,  is  explained  by  a  sharp  jump  of  tem¬ 
perature  from  °  =  6o°  to  r  =  65°  (see  Supplement  IV)  under  the  influence 
of  the  Gulf  Stream. 

If  the  conditional  equation  referring  to  the  65th  parallel  is  excluded 
from  the  compilation  of  the  normal  system,  the  value  of  the  mean  square 
error  will  decrease  while  by  will  remain  almost  unchanged  (see  solution 
(62)). 

The  comparatively  large  values  of  mean  square  errors  for  the  equations 
of  the  southern  hemisphere  in  the  winter  half-year  are  explained  by  a  non¬ 
uni  form  distribution  of  lands  and  oceans.  There  is  71  percent  of  land 
along  parallel  70°S,  and  not  more  than  9  percent  between  65°S  -  35°S.  If 
the  conditional  equation  pertaining  to  the  latitude  70°  is  excluded  from 
the  compilation  of  the  normal  system,  than  we  obtain  the  following  solu¬ 
tions 

(65) 


iiu,  r=  -5vG0  GW29C  U'4 
a'-—  :>.20-,-G.0027I'S:w  -• 

i* ’  •  v  •• 

u t  — G.S 1  1  v...i } 


As  it  appears  from  (59)  and  (6l),  the  coefficient  of  proportionality 
b8  varies  in  the  northern  hemisphere  from  0.0065  to  O.OOflU  and  in  the  south¬ 
ern  hemisphere  from  0.0071  to  0.0089.  It  Is  larger  for  land  and  smaller 
for  oceans.  There  is  almost  no  difference  between  values  ba  for  the  north¬ 
ern  and  southern  hemispheres.  On  the  average  the  value  of  the  coefficient 
bg  ■  0.0075  is  very  close  to  the  value  0.007,  determined  theoretically  by 
M.  Milunkovitch. 
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In  the  case  of  the  winter  half-year  the  value  of  the  coefficient  h* 
varies  from  0.003U  to  0.00^2  in  the  northern  hemisphere  and  from  0.0028  to 
0.0037  in  the  southern.  Here  there  is  also  no  difference  in  bv  for  the 
northern  and  southern  hemisphere;  in  the  winter  the  influence  of  the  dis¬ 
tribution  of  land  and  oceans  exerts  a  lesser  influence.  The  mean  value 
of  the  coefficient  of  proportionality  bw  is  0.0037*  i.e. ,  half  as  much 
as  for  the  fummer  half-year. 

In  the  course  of  three  million  years  of  the  geological  past  it  is 
possible  to  notice  an  increase  or  decrease  of  temperature  by  5  -  6°  in 
the  summer  and  2  -  3°  in  the  winter. 
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